Atomic Bose and Anderson glasses in optical lattices 
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An ultra cold atomic Bose gas in an optical lattice is shown to provide an ideal system for the 
controlled analysis of disordered Bose lattice gases. This goal may be easily achieved under the 
current experimental conditions, by introducing a pseudo-random potential created by a second 
additional lattice or, alternatively, by placing a speckle pattern on the main lattice. We show 
that for a non commensurable filling factor, in the strong interaction limit, a controlled growing 
of the disorder drives a dynamical transition from superfluid to Bose-glass phase. Similarly, in 
the weak interaction limit, a dynamical transition from superfluid to Anderson-glass phase may be 
observed. In both regimes, we show that even very low-intensity disorder-inducing lasers cause large 
modifications of the superfluid fraction of the system. 
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At zero temperature, the bosonic lattice gases may 
undergo a quantum phase transition Q from a super- 
fluid (SF) to an insulator 0- In the absence of disor- 
der there exist Mott insulator (MI) states, characterized 
by a fixed (integer) number of bosons per lattice site, a 
gap in the excitation spectrum, and vanishing superfluid 
fraction and compressibility. In the presence of disor- 
der, an additional insulating phase, so-called Bose-glass 
phase (BG), may occur Q. This phase presents no super- 
fluid fraction, a gapless excitation spectrum, and a finite 
compressibility. The SF-BG transition has been recently 
experimentally studied in various physical systems 3|, 
and attracts a continuous theoretical interest 0, 0, El 3 • 
In particular, the possibility of a direct MI-SF transition 
in the presence of disorder remains a controversial issue. 

The nature of the disorder-induced insulator phases 
depends on the interplay between hopping, nonlinear in- 
teractions and disorder. In the strong-interaction regime, 
the cooperation of interactions and disorder leads to 
the appearance of a BG phase. For weak interactions, 
the disorder leads to an Anderson- type insulator, or 
Anderson-glass (AG) 0. In the latter case, contrary to 
the BG phase, the interactions tend to delocalize and 
therefore compete with the disorder. 

The detailed analysis of these properties demands an 
experimentally accessible system in which the disordered 
Bose lattice gases could be studied in a controlled way. 
One of the aims of this paper is to show that this goal 
can be relatively easily accomplished by using cold Bose 
gases in optical lattices, for which the development of 
cooling and trapping techniques allows a large degree of 
control. Recently, Greiner et al. H, following the theo- 
retical suggestion of Jaksch et al. 0, have observed the 
SF to MI quantum phase transition in an optical (non 
disordered) lattice loaded with 87 Rb atoms. 

In the case of an optical lattice a pseudo-random po- 
tential can be dynamically generated by growing on an al- 
ready existing (main) optical lattice a second (additional) 



one with a different wavelength [ij. Hoi Hl| | . Based on this 
idea, we study in this Letter the dynamical generation of 
the BG and AG phases in a many-body system under re- 
alistic conditions. In addition, we discuss the alternative 
generation of a truly random potential by superimpos- 
ing a speckle pattern. In the first part of this Letter, 
we consider the regime of strong interactions where the 
filling factor is non-integer, showing that the adiabatic 
turn-on of the disorder may lead to a dynamical transi- 
tion into the BG phase. In the last part of the Letter, 
we analyze the weak interaction regime, and in particu- 
lar the possibility to achieve a dynamical transition into 
the AG phase. For both regimes of interactions, the SF 
fraction 0] is calculated. We show that under realis- 
tic conditions even very low-intensity disorder-inducing 
lasers lead to a dramatic reduction of the SF fraction, 
indicating the appearance of BG, or AG phases. 

We consider an ultracold Bose gas in a 2D optical lat- 
tice. We analyze the case of 23 Na, but our results also 
apply to other species. We assume the atoms as tightly 
confined in the transversal (z) direction by a harmonic 
trap of frequency lu z /2tt = 6kHz, so that the wavefunc- 
tion in z remains the Gaussian ground state. No addi- 
tional harmonic confinement is assumed in the xy plane. 

The optical lattice is formed by the main laser beams, 
whereas additional lasers are responsible for the introduc- 
tion of (quenched) pseudo-disorder. The optical potential 
has thus the following form: 



V{v) = Vi{v) + V r {v), 



(1) 



where r = [x,y], the main lattice is V;(r) = 
Vo(cos 2 (/ca;)-|-cos 2 (fcj/)), and the secondary one isV r (r) = 
yi(cos 2 (kir) + cos 2 (k2i")). The main (additional) beams 
intensity and/or detuning controls the value of Vb {V\). 
In the following we assume k\ = &2 7^ k, and the di- 
rections ki oc [—0.5,1] and k 2 cx [—1,0.5]. The pseudo- 
randomness of the potential is determined by q — k±/k, 
i.e. the ratio between the wavelengths of the main (A = 



2 



2ir/k) and the additional lattices. For commonly used 
NdYag and TiSapphire lasers q = 1064/795 = 1.338. We 
note that a pseudo random lattice can also be achieved 
by splitting off part of the main laser beams and creat- 
ing the additional incommensurate lattice by interfering 
these light beams at an angle. 

We assume that the energies involved in the system 
are much smaller than the energy separation between the 
first and the second band of the lattice, and consequently 
we can reduce our analysis to the first band. In that case, 
the physics of the atomic lattice Bose gas is governed by 
the Bose-Hubbard Hamiltonian |9j: 



H = - 



E 

<i,j> 



Jij a i a,j 



U< 



^^THK-lJ+^Wini. (2) 



Jq + 5 are the tunneling (hopping) coeffi- 

They slightly dif- 
by a correction of the form: SJij — 



Jij — 

cients between nearest neighbors 
fer from Jo 

— J d 3 rw*(r — ri)V r (r)w(r — rj), where w(r — rj) are the 
Wannier functions for the lowest energy band. Clearly, 
the contributions of the additional beams to the tunnel- 
ing, vanish on "average", < SJij 0. In the following 
we consider Vq = 25Er and V\ ~ 0.05.Er, where Er is 
the photon recoil energy. In this case we have checked 
numerically that \5Jij/ Jq\ < 0.1%, hence the major con- 
tribution to J.^ comes from J , and the model reduces to 
the ordinary Bose-Hubbard one with constant tunneling 
0,0,0,0- In the Hamiltonian |(2J, Uq oc a is the cou- 
pling constant [j| for the interparticle interactions (a is 
the scattering length), and 



tf, 



= J d 3 rw*(r 



ri)V r (r)w(r - r ; ) 



(3) 



are the pseudo-random on-site energies, which, as dis- 
cussed below, may introduce significant effects even for 
very small V\/Vq. 

The disorder-induced phases are characterized by a 
vanishing SF fraction, which is determined studying the 
system sensitivity to changes of boundary conditions. To 
this aim, we employ the boost method [13J, resulting in 



J ij 6 
(^2, Hi 



"*i . The angles ipij are 
and j — (xj,yj), then 



substitution in J2J: Jij 
defined as follows: if % = 
for yi = yj, <fij = sign(a;i - Xj)<p/M, and else ^ = 0, 
where M is the lattice size in the x-direction. Physi- 
cally, this choice of tpij corresponds to a constant current 
per lattice, proportional to (p, in the positive x direction. 
The SF fraction is then obtained from the corresponding 
ground-state energy E(ip) as 



M 2 E(<p) - E(0) 



N 



J f 2 



(4) 



where N is the number (mean) of atoms. Similar ex- 
pression can be easily derived for a ID case. In the fol- 
lowing, we denote by T the ratio of N to the number 



of lattice sites. The ground state is obtained from the 
minimization of {^mf\H — fiN\^ mf) , where \i is the 
chemical potential. We employ the Gutz wilier ansatz 
I*mf) = IliY,n \ n )ii wh ere fi l) are the ampli- 
tudes of having n atoms at a i-th lattice site. We have 
observed, that it is numerically easier to calculate the 
ground state using static methods for a non-disordered 
and non- twisted case (Vi = 0, <p = 0), and then adia- 
batically evolve such a state, first to ip ^ 0, and then 
for a constant <p to V\ ^ (first evolving V\ and then ip 
should give the same result). The evolution is performed 
by means of the dynamical Gutzwiller approach [l3 | , in 
which we solve the equations 



Ur 



n(n — 1) + nWi 



W + 



(5) 



where ^ = - E< lJ > Jif itpii {^MF\aj\$MF)- 

Another useful quantity characterizing the state of a 
Bose lattice gas is the condensate fraction, defined as 
the highest eigenvalue of the one particle density matrix, 
Pij = i^M f\ qIqi |^mf) i divided by the number of par- 
ticles 14]. This quantity is important for experiments, 
since it determines the phase coherence, and thus the 
contrast in interference measurements 8j. 
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FIG. 1: Condensate fraction (solid line) during the dynamical 
SF to BG transition induced by a superimposed additional 
optical lattice, for 23 Na atoms placed on the main 40 x 40 
lattice with A = 1064nm, q = 1.338, T = 0.75, V /E R = 25, 
and Vi(t)/E R = 0.059t/T where T = 0.42s is the total time 
of evolution. The inset shows the first preparation step (see 
text), for which Vo(t) = Vsf + {Vmi - V SF )t/T (Vi(t) = 0), 
where V S f/E r = 7, Vmi/Er = 25, and T = 0.02s. The 
dashed line presents the SF fraction. 



In order to study the dynamical transition into the 
BG phase, it is convenient to prepare the system in the 
SF phase in the presence of dominating interactions and 
weak tunneling, when T is non-integer. Therefore, we 
consider an initial system deeply in the SF regime, with 
almost 100% of SF fraction and T = 0.75. Then, the 
intensity of the main laser is adiabatically increased in 
20ms, obtaining a very large Uq/Jq w 70. Both, the 
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condensate and the SF fraction decrease during this pro- 
cess down to approximately 30%, being non-zero only 
due to non-integer value of T (see inset of Fig. QJ. As 
a next step, the disorder is turned-on adiabatically in 
about 0.5s, by switching on the additional laser beams. 
The condensate and also the SF fraction decrease dy- 
namically during this process (Fig. [I}. Ultimately, the 
condensate fraction does not tend to 0, but to a very 
small value of about 2%, due to the finite size of the 
systems and the approximate character of the Gutzwillcr 
approach 0,0]. In contrast, the SF fraction tends to 
zero faster. Note that the superfluidity is rapidly lost, 
although at any time the additional lattice is very much 
weaker than the main one. This fact could at first glance 
seem surprising, but it is due to the small values of Jo and 
Uq (10~ 3 Er and 7x 10~ 2 Er, respectively). Thus even for 
Vi being of the order of few percent of Er (V\/Vq <C 1), 
the value of |Wj|/E/o ~ 1. Together with the low value of 
Jo/Uo, this explains why the system enters the BG phase 
for such a weak additional lattice. 

We stress at this point that our calculations do not 
include an additional inhomogeneous trapping in the xy 
plane, which, for the low Jo we consider, may result in 
the formation of MI domains ^3] if the filling factor at 
the trap center is larger than 1. For sufficiently shallow 
traps, the central BG region should dominate the physics 
for a finite disorder. For filling factors lower than 1, the 
systems is expected to be fully in the BG phase. 
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FIG. 2: Condensate fraction during the dynamical SF to 
BG transition induced by a superimposed speckle pattern, 
for 23 Na atoms placed on the main 40 x 40 lattice with 
A = 1064nm, T = 1.5, Vq/Er = 25, and V 2 /E R = 0.048t/T, 
with T — 0.42s. The solid (dashed) line refers to speckles 
of average size F = 0.34A (r = 1.37A). For comparison, the 
dotted line shows the same transition in the quasi-disordered 
case (parameters as in Fig. Q but T = 1.5). We expect that 
similarly as in Fig. 0the decrease of the SF fraction is faster 
than that of the condensate one. 

It is interesting to compare results obtained with the 
quasi-disordered perturbation induced by the additional 
lattice, and those achieved using a purely random optical 
potential coming from a speckle pattern. We generated 
speckle pattern in a way described in |l8j and checked 
that it gives correctly both the autocorrelation function 



and the probability distribution. The speckle pattern 
induces a potential V s (f ) , which is characterized by its 
mean value U 2 = (V s (f)), and by the average speckle size 
r. As above, the disorder-induced corrections to the tun- 
neling are minor, and the most significant contribution 
of the speckle potential appears in the form of Wt coef- 
ficients |J3J, defined in this case by means of the V s (f) 
potential. Not surprisingly, values of Vi comparable to 
those previously considered for Vi, lead to a transition 
from the SF to Bose glass phase. We have investigated 
this issue for T = 0.34, 1.37, 2.75A. The obtained re- 
sults are similar to the ones resulting from quasi-random 
perturbation generated by additional beams (see Fig.|2). 
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FIG. 3: Transition from the SF to AG phase for 3 Na atoms 
placed on a ID main lattice of 86 sites, with A = 1064nm 
and Vo/Er — 16. Upper plot: superfluid fraction vs. ampli- 
tude of the disordered potential. The thick curve indicates 
the speckle-induced transition (r = 1.34A), whereas the nor- 
mal line refers to the transition generated by a superimposed 
lattice with ki/k = 1.338, k 2 /k = 1.396 (see text). Middle 
plot: dependence of V r (x) (normal line) and V a (x) (thick line) 
for Vi/Er = 0.004 and V 2 /E R = 0.001. Lower plot: occupa- 
tion numbers of the ground state of the Hamiltonian JSJ in 
the presence of speckles (filled circles) and additional beams 
(hollow circles), for Vi/Er = 0.004 and V 2 /E R = 0.001. 



In the final part of this Letter, we discuss the non- 
interacting regime, where the adiabatic turn-on of the 
disorder can lead to a dynamic transition into the AG 
phase. The weakly interacting regime, and even the 
non-interacting one, could be achieved by reducing the 
s-wave scattering length by means of Feshbach reso- 
nances [l^. For simplicity we consider a ID opti- 
cal lattice case where the perturbation is either gener- 
ated by two non-commensurate standing waves, V r {x) — 
Ui(cos 2 (A:ix) + cos 2 (fc 2 x)) with k\ 7^ fe, or by a laser 
speckle potential V s (x) characterized by a mean value 
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Vi- The Hamiltonian JSJ reduces then to the form 

H = ~ ^ Jij a i a 3 + X! W i n i> ( 6 ) 
<i,j> i 

which becomes the famous Anderson's hopping model 
[20| . if the Wi coefficients are randomly distributed. Sim- 
ilarly as above Jij are almost unchanged by the pres- 
ence of disorder. The single-particle character of the 
non-interacting problem greatly simplifies the calcula- 
tions and allows for an exact treatment. In the absence 
of interactions, 100% of atoms condense at zero temper- 
ature, each in a single particle state |^) = ^2 n c n ajjvac). 
We have analyzed the dynamical evolution of the ampli- 
tudes c„ during the turn-on of the disorder as well as 
their ground state distributions — see Fig. |3J 

Due to the finite size of the system, which implies dis- 
creteness of the excitation spectrum, the application of 
the boost method in the absence of disorder results in the 
100% SF fraction. When the typical energy of disorder 
becomes comparable to the energy separation between 
the energy levels in the absence of disorder, the SF frac- 
tion is expected to vanish (see Fig. We clearly observe 
a transition from a delocalized SF state to an Ander- 
son localized one, where the occupation probabilities of 
neighboring sites decrease exponentially with the distance 
from the localization centers [2(J- From the upper plot 
in Fig. |21 we conclude that in order to drive the system 
into the AG phase, smaller intensities are needed for the 
speckle pattern than for the superimposed incommensu- 
rable lattice, since as expected the speckle pattern results 
in a "more disordered" distribution of the Wi coefficients 
|2l| (see middle plot in Fig. |3J|. It is especially worth 
to stress that in the case of a superimposed incommen- 
surable lattice the periodicity of the lattice reflects itself 
in periodically Anderson-localized domains, but the local- 
ization phenomena is still present, as shown in Fig. [3] We 
observed the same quantitative phenomena also for ratios 
ki/k and k2/k different than the ones listed in the caption 
of Fig. We have performed also time-dependent simu- 
lations based on Eq. © to determine the time scale for 
the adiabatic transition into the AG phase. For a main 
lattice wavelength A = 1064nm the adiabatic evolution 
lasts few seconds. This may be shortened using higher 
frequency lattice beams since evolution time scales as A 2 . 

The experimental observation of the SF to BG (or to 
AG) transition should be relatively easy to accomplish in 
a setup as that of Ref. ||- The BG and AG phases can 
be detected by observing the interference pattern after 
removing the lattice for different intensities of the super- 
imposed laser beams. The insulator character of the BG 
and AG phases will be revealed by the disappearance of 
the interference fringes. The phases can be additionally 
characterized by measuring their gapless excitation spec- 
trum, in an experiment similar to that of Ref. 

Summarizing, we have proposed an experimentally fea- 



sible and relatively simple method of creating a system 
whose physics is governed by the disordered Bose Hub- 
bard model. We have shown how the onset of small per- 
turbation of the lattice potential may result in a dynami- 
cal transition from a supcrfluid regime into Bose-glass, or 
Anderson-localized phases. This transition occurs within 
experimentally feasible time scales and can be easily con- 
trolled, allowing for a detailed analysis of disorder in- 
duced transitions. Our proposal stimulates in this sense 
new interesting experimental possibilities including stud- 
ies of Anderson localization in 2D systems (which is still 
a controversial topic), and the investigation of the SF to 
MI transition in the presence of disorder. 
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